Production of (j) Mesons in Au + Au Collisions at 11.7AGeV/c 
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We report on a measurement of 0-meson production in Au + Au collisions at a beam momentum 
of 11.7AGeV/c by Experiment 917 at the AGS. The measurement covers the midrapidity region 
1.2 < y < 1.6. Transverse-mass spectra and rapidity distributions are presented as functions of 
centrality characterized by the number of participant projectile nucleons. The yield of 0's per 
participant projectile nucleon increases strongly in central collisions in a manner similar to that 
observed for kaons. 

PACS numbers: 25.75.-q, 25.75.Dw, 13.85.Ni 



I. INTRODUCTION 

The production of mesons in relativistic heavy-ion 
collisions has been an important subject of study at 
the AGS, the SPS, and RHIC. The production of the 
(j) meson, the lightest bound state of strange quarks 
(ss), is suppressed in ordinary hadronic interactions be- 
cause of the Okubo-Zweig-Iizuka (OZI) rule [Q. It has 
been proposed that in a quark-gluon plasma (QGP) sce- 
nario, strange quarks could be rapidly and abundantly 
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produced via giuon interactions Thus, ip mesons 

could be created in a non-conventional way via strange 
quark coalescence, bypassing the OZI rule. A strong en- 
hancement in 0-meson production would serve as one of 
the strangeness-enhancement signatures for QGP forma- 
tion |. 

The production of (f> mesons has been measured in 
Si + Au coUisions at 14.6AGeV/c by E859 at the AGS 
within a rapidity interval of 1.2 < y < 2.0 Q. The ra- 
tio of the total cj) yield to the yield has been found 
to be 10% for the uppermost 7% of the charged-particle 
multiplicity distribution. NA49 at the SPS has measured 
(/)-meson production in p + p, p + Pb, and Pb -I- Pb col- 
lisions at a beam energy (i^bcam) of 158A GeV within a 
rapidity range of 3.0 < y < 3.8. An enhancement in the 
ratio of to pion total yields by a factor of 3.0 ± 0.7 has 
been observed in central Pb -I- Pb collisions relative to in 
p + p interactions §. Also at the SPS, NA50 has mea- 
sured the ratio ip/ip + as a function of centrality in 
Pb -|- Pb collisions by fitting the invariant-mass spectrum 
of muon pairs with < y ~ yNN < 1 and 1.5 < mf < 3 
GeV/c^ Both this ratio and the absolute yield of 0's 
per participant nucleon have been found to increase sig- 
nificantly with the number of participant nucleons. The 
STAR experiment has reported on a measurement of the 
(f) yield at midrapidity in Au -I- Au collisions at center-of- 
mass energy y^SNN = 130 GeV. An increase in the (p/h^ 
ratio with ^/snn has been observed . 

Enhancements in the yields of particles with open 
strangeness have also been observed. At the AGS, 
K^/tt^ ratios have been seen to increase to ^20% in 
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central Si + Au and Au + Au collisions, from 4-8% in 
p + p collisions [|[ ^, |l^, and the A/p ratio has been 
found to increase strongly with centrality in Au + Au 
collisions At the SPS, the WA97 collaboration 

has observed enhanced production of K's, A's, S's, and 
ri's in heavy-ion collisions relative to in p + p or p + A 
collisions 13 1. Measurements from the PHENIX exper- 
iment at center-of-mass energy ^/snn = 130 GeV have 
shown the kaon yield to be more strongly dependent on 
centrality than the pion yield at midrapidity [ij] . 

The reason for the strangeness enhancement in heavy- 
ion collisions is not completely understood. The rescat- 
tering of hadrons and the conversion of the excitation en- 
ergy of secondary resonances into strange particles might 
give rise to strangeness enhancement in a purely hadronic 
picture (l|, |l|, |l|. Stran geness enhancement can also 
be interpreted as a reduction in canonical strangeness 
suppression from p +_p to A + A reactions in the context 
of thermal models [[l8| . Strangeness enhancement as a 
function of the number of "grey" protons (a centrality 
index) has also been seen in p -|- A collisions, in which a 
QGP phase is unlikely to contribute [ p^ . 

Under the conventional hadronic interactions, there are 
mainly three scenarios proposed for the production of 
mesons in nuclear collisions: 

• Parton fusion of strange sea quarks or knock- 
out of ss pairs pl| from the primary collisions of 
projectile and target nucleons. 

• Secondary baryon-baryon interactions BB 
(j)NN, meson-baryon interactions {tt,p)B 

4>B ll^l and meson- meson interactions np (j) , 
where B ^N,A,N*. 

• Secondary kaon-hyperon interactions KY — > (j)N 
and kaon-antikaon scattering KK (f>p in the 
event of the restoration of chiral symmetry in the 
hot and dense nuclear fireball |Q. 

With inclusion of secondary meson-baryon and meson- 
meson interactions, RQMD (Relativistic Quantum 
Molecular Dynamics) is able to qualitatively describe 
the increase of kaon yields with centrality measured by 
E866 1^ , though quantitative differences do exist |^ . 

The different scenarios described above imply different 
relations between the scaling of </> production and that 
of other hadrons, such as kaons and pions, with collision 
centrality and energy. Hence, a systematic measurement 
of (p production in different collisional systems may help 
to quantify the increase in the overall strangeness pro- 
duction, and, in combination with other measurements 
of strange and non-strange hadron production, to differ- 
entiate between possible mechanisms contributing to the 
strangeness enhancement [ p7| . 

In this paper we report on a measurement of the 
yield around midrapidity in Au -I- Au collisions at the 
Alternating Gradient Synchrotron (AGS) at Brookhaven 
National Laboratory (BNL), and compare it to the yields 
of pions and kaons. 



II. EXPERIMENTAL DETAILS 

Experiment 917 took data on Au+ Au reactions at pro- 
jectile momenta of 6.8, 8.9 and 11.7AGeV/c (correspond- 
ing to center-of-mass energies ^snn = 3.83, 4.31, and 
4.87 GeV, respectively) in the fall of 1996. The identi- 
fied particles include particles containing strange quarks, 
such as K^, K~ , 0, A, and A, and non-strange particles 
such as TT, p, and p ||l^, |l^ ^ |2^, The experimental 
apparatus consisted of a movable magnetic spectrometer 
for tracking and particle identification, and beamline de- 
tector arrays for global event characterization. When the 
data reported in this work were collected, the beam mo- 
mentum was 11.7^4 GeV/c, an Au target with areal den- 
sity 1961 mg/cm^ (corresponding to approximately 4% 
of an interaction length for an Au projectile) was used, 
and the spectrometer angle was set to either 14° or 19° 
from the beam axis. The momentum resolution of the 
spectrometer, Sp/p, was about 1% for particles with mo- 
mentum greater than 1 GeV/c and increased at lower mo- 
mentum, up to 2% at p = 0.6 GeV/c, due to the effect of 
multiple scattering. The kaon momenta in reconstructed 
(j) events lay mostly in the range of 1.0 < p < 2.0 GeV/c. 
Event centrality was characterized by the energy of the 
beam spectators after the interaction as measured in the 
zero-degree calorimeter (ZCAL), positioned downstream 
of the target on the beam axis. More details on the de- 
tector systems are given in Refs. ^ 

The data presented here were collected using a 
two-level online trigger: a minimum-bias spectrometer- 
activity trigger (LVLl) followed by a particle- 
identification trigger (LVL2) which required two 
charged kaons of either sign or one p in the spectrometer 
acceptance. The hardware LVL2 trigger looped over the 
combinations of drift chamber and TOE hits, and formed 
combinations consistent with tracks corresponding to 
particles of given momentum, charge, and particle type 
using a look-up table |Q. The LVL2 trigger increased 
the live time of the data-acquisition system by essentially 
the ratio of the rate of vetoable LVLl triggers to the 
rate of vetoable LVLl triggers that were not vetoed. 
This ratio was known as the LVL2 rejection factor, and 
could be estimated online. In Au + Au collisions with 
a magnetic field of 4 KG and the spectrometer at 14° 
and 19°, typical rejection factors for a 2K/p trigger were 
4.8 and 8.5, respectively. By examining data that was 
taken with the LVLl trigger and recording the decision 
of the LVL2 trigger, the inefficiency of the latter was 
monitored and found to be less than 1% for events fully 
inside the acceptance. Most of the LVL2-triggered data 
was background since the trigger was optimized to reject 
only events that were clearly not of the correct type in 
order to keep its efficiency near 100%. 

A time-of-fiight (TOE) system with a typical resolu- 
tion of 130 ps served to identify pions and kaons up to a 
momentum of 1.75 GeV/c with 3 standard deviations of 
TOE resolution. Above this momentum, the 3a contours 
in the (p, TOE) plane began to overlap. Within this re- 
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TABLE I: Centrality bins used in the analysis. The cuts on 
zero-degree energy, -Ezcal, that define each bin are fisted 
(in GeV), together with the corresponding fraction of the 
total cross section for Au + Au collisions (in percent), the 
mean number of projectile participants for the bin, (Npp), 
and the estimated value of the mean number of binary col- 
lisions, (Ncou), by Glauber model and the mean number of 
projectile participants estimated from -Ezcal, (Afp^^'^). The 
total beam kinetic energy in the collisions is about 2123 GeV. 



bin 


SzCAL 


% CTtot 




(Afcoll) 




1 


0-280 


0-5 


170 ± 10 


770 ± 72 


179 ± 1.0 


2 


280-560 


5-12 


138 ± 12 


585 ± 74 


157± 1.5 


3 


560-960 


12-23 


101 ± 14 


393 ± 74 


126 ± 2.0 


4 


960-1440 


23-39 


62± 14 


207 ± 60 


84 ± 2.6 


5 


>1440 


39-77 


20± 12 


50 ±38 


34 ±3.1 



gion, particles were identified as kaons only if their TOF 
was inside the kaon region and outside of the pion region, 
since 10-30% of the tracks in the overlap region were from 
pions. A specific correction for the kaon- identification 
inefficiency from the exclusion of the pion bands in this 
region of high momenta was implemented in the analysis 
(see Sec. [II B). Additionally, kaons with momenta less 
than 0.5 GeV/c were rejected, in order to avoid the need 
for large acceptance corrections. 

The index of event centrality was obtained from the 
distribution of the energy deposited in the zero-degree 
calorimeter, -Ezcal, obtained using a minimum-bias 
interaction trigger based on the measurement with a 
Cerenkov counter of the projectile charge after passage 
through the target. Empty-target runs were periodically 
taken in order to subtract the background from the ZCAL 
energy distribution. Assuming a monotonic relation be- 
tween the energy deposition in the calorimeter and the 
event centrality such that the most central events corre- 
spond to the smallest energy deposition, we divided the 
data into five centrality bins. The centrality bins are ex- 
pressed as a fraction of the total cross section for Au±Au 
collisions, 6.8 b as evaluated from the parameterizations 
of Ref . . The cross section for Au ± Au collisions as 
measured with the minimum-bias interaction trigger is 
approximately 5.2 b, or 77% of the total cross section. 

We have attempted to estimate (Npp), the mean num- 
ber of projectile participants, and (A^cou)) the mean 
number of binary collisions, according to the Glauber 
model using impact parameter cuts correspondini 
to the centrality range for each bin as listed in Table 
and an assumed p + p inelastic cross section of 30 mb. In 
Experiments 802, 856 and 866 the number of projec- 
tile participants has been conventionally estimated from 

(1 - E'zCAL/^^bcam), whcrC 



the relation iVp^'^'" = Ap^oj 



iproj 



is the atomic mass of the projectile, and -Ebeam 
is the total kinetic energy of the nuclei in the beam be- 
fore interaction, which was 2123 GeV in the present case. 
For the purpose of comparison, we have also calculated 
Npp^^ and its uncertainty with a nominal ZCAL energy 
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FIG. 1: The experimental acceptance for K'^ K~ pairs in the 
space of transverse mass {mt) vs. rapidity (y) for the 19° and 
14° spectrometer angle settings. 



resolution [aE/E) of 3.6% |Q. All relevant parameters 
for each centrality bin are listed in Table |. To facilitate 
comparison with other experiments we will use {Npp) as 
calculated using the Glauber model for the rest of this 
paper. 



III. EXPERIMENTAL RESULTS 

A total of 250 million LVL2-triggered events were an- 
alyzed. When the two spectrometer angle settings were 
combined, the overall acceptance for K^K^ pairs cov- 
ered the region of rapidity 1.2 < y < 1.6 and trans- 



verse mass 1.0 GeV/c^ < 



mt 



{Ptic) 



211/2 



< 



2.2 GeV/c^, as shown in Fig. |. 

The (j) mesons were reconstructed by forming the 
invariant-mass (rriinv) distribution of identified K^K~ 
pairs and subtracting the combinatorial background, 
which was obtained by the event-mixing method [ p7[ . 
Two kinds of distributions were formed: "same-event" 
distributions, in which the X+'s and if ^'s were selected 
from the same event, and "mixed-event" distributions, 
in which the individual particles were chosen from differ- 
ent events in the same centrality class to represent the 
uncorrelated background. 



A. Invariant- mass distribution for K pairs 

First, we examined the minimum- bias invariant-mass 
distribution of kaon pairs by fitting the data to a relativis- 
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tic Breit-Wigner distribution (RBW) convoluted with 
a Gaussian representing the experimental mass resolu- 
tion 1^, The shape of the background (BG) was ob- 
tained from the mixed-event distribution while the nor- 
malization was left as a free parameter in the fit. The 
parameterization is as follows: 



J /\rsame 

dm 



rm2 exp 

a / RBW(m') 

J mi 

J A r mixed 



1 ( m—rn \2 
" 2 V cr^ ) 



dm! 



dm 



where 



RBW(m) 



mmor(m) 



r(m) = 2rn 



(m? — mf^Y + (mor(m))2 ' 



(9/90)2 + 1' 



90 



\J mQ /4 — m 



2 

AT' 



(1) 

(2) 
(3) 

(4) 
(5) 



and where the limits of integration were mi = 
0.989 GeV/c^ and = 1.252 GeV/c^. There were five 
free parameters in all: the relative normalizations a and 
6, the peak mass and width parameters toq and To, and 
the experimental mass resolution (t,„. 

For the minimum-bias data, the same-event invariant- 
mass distribution of K'^ K~ pairs and the background- 
subtracted signal distribution corresponding to the first 
term in Eq. (|l|) are shown together with the fits in 
Fig. |. The values toq = 1018.99 ± 0.36 MeV/c^, 
Fo = 6.14 ±2.59 MeV/c^, and cr,„ = 2.43 ±1.11 MeV/c^ 
were obtained with xV^of = 145/167 = 0.87. The fit 
results for the peak position (mo) and width (Fq) of the 
(j) signal are in agreement with the world-average values. 
The large uncertainty on Fq arises mostly because of a 
high degree of correlation in the fit between the values 
for the parameters Fg and h (correlation coefficient = 
0.71). The value obtained for the experimental mass res- 
olution {am) is also consistent with our estimated value 
of 2.0 MeV/c^, from Monte Carlo studies and from the 
width of the peak from A decay in the p7r~ invariant- 
mass distribution (the known contribution from the 
multiple scattering of the kaons in the target in the case 
of the present measurement was accounted for by addi- 
tional smearing). It is noted that the result of the fit 
lies systematically below the data in the invariant-mass 
region 10-20 MeV/c^ below the peak. We do not have 
a clear understanding of this effect. 



B. Signal counting and extraction of differential 
yields 

Monte Carlo studies confirm that the experimen- 
tal invariant-mass resolution {<Jm) remains constant 
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FIG. 2: (Color Online) The invariant-mass {rriinv) distribu- 
tion of K~ pairs for minimum-bias events, superimposed 
witfi a fit consisting of a relativistic Breit-Wigner distribution 
plus the background distribution from a mixed-event tech- 
nique, as described in the text. In the inset, the rf) signal to- 
gether with fit results from the first term in Eq. (nl) is shown. 



(^2 MeV/c^) across the kinematic region of accep- 
tance Q. Given the stable mass resolution and the 
limited statistics for the division of data into bins in cen- 
trality and phase space, the transverse-mass spectra were 
obtained by counting events in a defined signal window 
and estimating the number of background events within 
this window, rather than by integration of the result of 
the fit withEq. (0). The details were as follows. The data 
were divided into subsets for each of the two spectrom- 
eter settings and five centrality classes. Same-event and 
mixed-event pairs for each subset were binned in rapidity 
and transverse mass (Ay = 0.2 and Am* = 0.2 GeV/c^), 
and invariant-mass distributions were obtained for each 
bin. No corrections were applied to the invariant-mass 
distributions at this stage — signal extraction was per- 
formed on the spectra of raw counts. In each bin, the 
total number of signal and background counts {S + B) 
was taken to be the sum of the counts within a signal 
window defined as 1.0040 < mi„„ < 1.0355 GeV/c^. In 
units of the natural cf) width, this interval corresponds to 
approximately ±3.7Fo about the expected peak position. 
We then estimated the number of background counts {B) 
inside the signal window as the number of counts in the 
normalized distribution of mixed-event pairs in this same 
interval. For this purpose, the mixed-event distributions 
were normalized to have the same area outside of the 
signal window as the same-event distributions. The final 
number of signal counts {S) was obtained by subtract- 
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ing the background estimated in this way from the total 
number of events in the signal window. Using this proce- 
dure, we were able to obtain estimates of S for a total of 
88 distributions corresponding to different bins in phase 
space and centrality for each of the two angle settings. 

We use the statistic obtained from the Poisson log- 
likelihood treatment described in Rcf. js^ as an index of 
the degree of consistency between the same- and mixed- 
event invariant-mass distributions in the mass region 
used for normalization. We typically found x^/dof « 0.8 
(for most fits, the number of degrees of freedom — one mi- 
nus the number of populated bins outside of the signal 
window in the mixed-event distributions — was 149). In 
the worst case, we found x^/dof = 177/149. 

The statistical error on S was calculated as cr|. = 
{S + B) + a'^. Here, (5 -I- B) is the contribution from 
Poisson fluctuations in the number of total counts in the 
signal window, and ctb was taken to be B/^/B^, with 
Bout the number of counts in the spectra outside the sig- 
nal window (where same- and mixed-event spectra are 
normalized to the same integral). As an index of the 
statistics of the measurement, the values obtained for S 
range from a few counts for bins in phase space near the 
limits of the spectrometer acceptance, to approximately 
100 counts for bins well within the acceptance. Typical 
values of the signal-to-total ratio, S/{S + B), were 0.05- 
0.30; 74 of the 88 measurements had values of S/{S + B) 
within this interval. The median statistical error for the 
88 measurements of S was 72%; 57 of the 88 measure- 
ments had statistical errors of less than 100%. 

Systematic errors on S may result if the signal window 
does not completely include the signal peak. In this case, 
not only are the tails of the signal peak excluded from 
the signal count — the equal-area normalization technique 
will also cause the baseline to be overestimated. System- 
atic errors on S may also arise if the mixed-event dis- 
tributions do not accurately represent the background in 
the same-event distributions. Of particular concern is the 
excess of cou nts in the mass interval 1.000-1.010 GeV/c^ 
noted in Sec. [II A| and visible in Fig. |[ In the same-event 
spectra for the individual centrality and phase-space bins, 
this excess is present to varying degrees. Since we are 
uncertain of the origin of this excess, its effects must be 
included in our estimates of the systematic error for the 
signal count in each bin. 

We assign the systematic error on S for each bin by 
varying the width of the signal window from ±2ro to 
±5ro. For window sizes smaller than ±2ro, the signal 
count drops precipitously. For window sizes larger than 
±5ro, there is very little constraint on the normalization 
from the low-mass portion of the distribution. In Fig. ^ 
we show the dependence of S on the width of the sig- 
nal window for minimum-bias distributions in six bins 
of y and mj. For each bin, values of 5' are plotted for 
nine different choices for the signal window, each sym- 
metric about the expected peak position with half-width 
in units of Tq given by the abscissa value. The central 
point in each plot, corresponding to a window width of 
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FIG. 3: Dependence of S on width of signal window for six 
bins in y and mt: la) 1.2 < y < 1.4, 1.5 < mt < 1.7 GeV/c^, 
2a) 1.2 <y < 1.4, 1.7 < rut < 1.9 GeV/c^ 3a) 1.4 < y < 1.6, 
1.7 < mt < 1.9 GeV/c^ lb) 1.2 < y < 1.4, 1.3 < mt < 
1.5 GeV/c^ 2b) 1.4 < y < 1.6, 1.3 < mt < 1.5 GeV/c^ and 
3b) 1.4 < t/ < 1.6, 1.5 <mt < 1.7 GeV/c^ For panels la, 2a, 
and 3a, the spectrometer was at 19°; for panels lb, 2b, and 
3b, the spectrometer was at 14°. Minimum-bias data (0-77% 
central) are shown in each case. The abscissa is the half- 
width of the signal window expressed in units of the (j> width, 
To = 4.26 MeV/c^. The errors plotted are statistical only. 
The lines represent the results of linear, least-squares fits to 
the points, using uniform weights, and are used to obtain the 
estimate of the systematic error on S for each bin. 



about ±3.7ro, provides the value of S used to obtain the 
results presented in this work. Some of the point-to-point 
fluctuation is statistical in nature and is presumably ac- 
counted for in the statistical errors associated with each 
value of 5. In most cases, however, there is a sugges- 
tion of an underlying systematic trend. To quantify this 
trend, we perform a linear, least-squares fit with uniform 
weights to the nine points for each bin. The results are 
shown as the straight lines in the figure. For each bin, 
we quote as the systematic error on S the RMS varia- 
tion of this line over the interval spanned by the first and 
last points, which is equal to the maximum vertical ex- 
tent divided by ^/l2. When applied to each centrality 
and phase-space bin used for the analysis, this procedure 
leads to systematic errors for S that vary significantly 
from bin to bin. For the 88 individual measurements of 
S, the median systematic error is 16%; 62 of the 88 mea- 
surements have systematic errors of less than 30%. In 
all cases, the statistical errors on S are larger than the 
systematic errors. 

The differential yields for each bin in phase space and 
centrality were then calculated from S as 



1 



1 



^ S {a){w) 

dydmt AyAmt /zcal^^int BR(mt 



(6) 



where /zcal is the fraction of minimum-bias interaction 
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triggers satisfying the centrality selection, A'int is the 
background-subtracted number of minimum-bias inter- 
action triggers collected with the data for normalization 
purposes, (a) is the geometrical acceptance correction, 
(w) is the overall correction factor for other experimen- 
tal effects, (mt) is the mean value of mt for all pairs in 
the bin, and BR{(f> -> K+R-) = 49.2%. 

Apart from the systematic error on S, the principal 
source of systematic uncertainty in the yields is from 
the uncertainty on /zcal- Due to radiation damage to 
the plastic-scintillator materials of the ZCAL, its signal 
decreased with time. A run-dependent calibration was 
therefore used, in which the energy of the fragmented- 
beam peak was recentered at the nominal beam kinetic 
energy (with corrections for energy loss in the target). 
These calibrations were performed periodically through- 
out the running period. The centrality bins were then 
defined via fixed cuts on the calibrated ZCAL energy. 
However, /zcal exhibited residual variation as a function 
of running time. The resulting contributions to the sys- 
tematic errors on the yield measurements are estimated 
to range from about 5% in the most central bin to 10% 
for the most peripheral bin. 

The geometrical acceptance, (a), and many of the con- 
tributions to the weight factor {w) , were calculated using 
a GEANT-based Monte Carlo simulation implementing 
a realistic detector configuration js^. The value of (a) 
used in Eq. (^) is the average value for </> — ^ K^K^ 
events in the phase-space bin of interest. The weight fac- 
tor {w) was calculated as follows for each centrality and 
phase-space bin. Corrections for kaon decays in flight, 
the single-particle tracking inefficiency, hadronic interac- 
tions, multiple scattering in the spectrometer, momen- 
tum cuts, and the kaon-identification inefficiency were 
evaluated using the Monte Carlo, separately for each 
track in each pair. An additional correction for the single- 
track tracking efficiency in the presence of background 
hits on the drift chambers was evaluated by inserting 
found tracks into random events. This occupancy cor- 
rection was then parameterized by the hit multiplicity in 
the chambers and applied for each track; the resulting 
correction was typically 20-40%. The individual correc- 
tion weights for each track and from each of the above 
sources were then multiplied, together with an additional 
correction for a two-track opening-angle cut (evaluated 
for each pair), to obtain the overall pair weight, w. The 
value of {w) in Eq. (^) is the average value of w for all 
pairs in the same-event distribution of interest. The sys- 
tematic errors associated with (a) and {w) are estimated 
to be no larger than 5%, arise mainly from the statistics 
of the samples used to obtain the weight factors, and do 
not significantly affect the results. 

The final systematic errors on the individual yield mea- 
surements are calculated as the quadrature sum of the 
fractional systematic errors on S (for each measurement), 
on /zcal (for each centrality bin), and on {a){w) (con- 
stant). In general, the dominant contribution is from the 
systematic uncertainty on S. For the 88 individual differ- 



ential yield measurements, the median overall systematic 
error is 20%; 60 of the 88 measurements have systematic 
errors of less than 30%, and 75 of the 88 have systematic 
errors of less than 50%. In all cases, the statistical er- 
rors are larger than the overall systematic errors; for 73 
of the 88 measurements, the statistical errors are more 
than twice as large as the overall systematic errors. 

Experimental results from both the 14° and 19° set- 
tings of the spectrometer were in all cases consistent 
within statistical errors, and were therefore combined for 
the presentation of the transverse-mass spectra, result- 
ing in the 53 individual differential yield measurements 
plotted in Fig. ^. 



Transverse-mass spectra and rapidity 
distributions 



The transverse-mass spectra for the five centrality bins 
are shown in two bins of rapidity, 1.2 < y < 1.4 and 
1.4 < y < 1.6, with an mt bin size of 0.2 GeV/c^, in 
Fig. ^. For the rrit bins where the acceptance from 14° 
and 19° data sets overlap, the individual measurements 
from each of the data sets agree to within statistical er- 
rors, and the weighted average of the two measurements 
is presented. The full error shown includes both statisti- 
cal and systematic contributions; the contribution from 
statistics alone is indicated by the cross bars. The to- 
tal errors are calculated as the quadrature sum of the 
statistical and systematic errors. 

An exponential parameterization with two parameters, 
rapidity density (dN/dy), and inverse slope (T), was used 
to fit the transverse-mass spectra: 



1 



dN/dy 



2'KC'^mt dmtdy 



27r(TTOoc2 + T2) 

(mt - mo)c2 



X exp 



T 



(7) 



This function gave a good fit to the data with 
X^/dof^l in all cases. The values obtained for the rapid- 
ity density (dN/dy) and inverse-slope parameters (T) are 
plotted as functions of rapidity and centrality in Fig. |^ 
and are listed in Tables || and III. The statistical and 
systematic contributions to the error were evaluated by 
performing the fits with and without the systematic con- 
tribution included in the errors. The systematic error due 
to uncertainties in the overall normalization is excluded 
from the systematic error estimate for T. 

As a function of centrality, the rapidity density shows 
a strong systematic increase, while the inverse slope in- 
creases more mildly. Within the rapidity range covered 
by the measurement, there seems to be no strong rapidity 
dependence for T. 

The fiducial yield is the sum of the values dN/dy x Ay 
for the two rapidity bins covering the interval 1.2 < y < 
1.6 and is plotted as a function of (Npp) (calculated from 
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FIG. 5: The inverse-slope parameters, T, and rapidity den- 
sity, dN/dy, of (j) mesons as a function of centrality. The open 
symbols show the reflections of the data points about midra- 
pidity (v — 1.6). The error bars have the same signiflcance as 
in Fig. 



FIG. 4: The invariant yield of 0's as a function of transverse 
mass (mt — mo) in five centrality bins labeled by a/atot- In 
each panel, the top set of points (circles) is for 1.2 < y < 
1.4, while the second set of points (triangles) is for 1.4 < 
y < 1.6 and is divided by ICQ for presentation. The total 
errors are calculated as the quadrature sum of the systematic 
and statistical errors, while the magnitudes of the latter are 
indicated by the cross bars. The lines are the results of the 
exponential fits described in the text. 



a Glauber model as explained in Sec. ||) in Fig. |^. When 
normalized to (iVpp), this quantity exhibits a steady in- 
crease with increasing centrality. This implies that the 
fiducial yield of (/<'s increases faster than linearly with 
(iVpp). It is noted that the increase with centrality is 
slightly more significant (the effect is on the order of 4% 
in the slope of a linear fit) if the fiducial yields are plot- 
ted against (N^^^^). Thus, we observe an enhancement 
in 0- meson production in central events. This is similar 
to what has been previously observed for kaon produc- 
tion @. In Pb -f Pb coUisions at the SPS, NA50 has also 
observed an increase in the yield of (/)'s per participant 



nucleon ((np)) as a function of ( 



] , though the NA50 



data show more evidence of saturation of this quantity 
in the most central collisions than do our results. 



IV. DISCUSSION 



In this section, we compare the measured inverse slopes 
(T) and rapidity densities {dN/dy) for the (j) to those 
of the other hadrons and explore their dependence on 
centrality and beam energy. 



TABLE II: Inverse-slope parameter, T, in units of MeV/c^, 
for <^'s in two bins of rapidity and five bins of centrality. Errors 
are statistical followed by systematic. 



Centrality 



1.2 <y < 1.4 



1.4 < 2/ < 1.6 



0-5% 

5-12% 

12-23% 

23-39% 

39-77% 



196 ± 19 ± 19 
244 ± 20 ± 22 
161 ± 18± 19 
154 ±15 ±21 
169 ± 17 ± 15 



260 ± 27 ± 25 
173 ± 28 ± 50 
215 ±28 ±25 
244 ± 31 ± 40 
189 ± 28 ± 90 



Transverse-mass spectrum of the 
to other species 



comparison 



In the most central bin (cr/atot = 0-5%), the inverse- 
slope parameter is determined to be 260±37 MeV/c^ 
around midrapidity. As motivated by hydrodynamic 
models |Q, the inverse-slope parameters for different 
particle species are expected to scale with species mass 
for particles that participate in the collective transverse 
fiow. It is interesting to see whether the inverse slope for 



TABLE III: Rapidity density, dN/dy, for (^'s in two rapidity 
bins, by centrality. Errors are statistical followed by system- 
atic. 



Centrality 



1.2 <y < 1.4 



1.4 < y < 1.6 



0-5% 0.312 ± 0.075 ± 0.056 

5-12% 0.203 ± 0.032 ± 0.018 

12-23% 0.124 ±0.040 ±0.020 

23-39% 0.098 ± 0.023 ± 0.018 

39-77% 0.017 ± 0.005 ± 0.004 



0.362 ± 0.065 ± 0.054 
0.316 ±0.122 ±0.153 
0.158 ±0.045 ±0.069 
0.076 ±0.017 ±0.022 
0.018 ±0.005 ±0.005 
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FIG. 6: The fiducial (f} yield (left) and fiducial yield per 
projectile participant, (Npp), as a function of (Npp) (right). 
The accepted rapidity region is 1.2 < 1.6. The error bars 
have the same significance as in Fig. 



(j) mesons fits into the systematic trend observed for the 
other hadrons. However, two known effects complicate 
the interpretation. 

The first is that the cross sections for cf) interactions in 
the collision medium are expected to be small. The OZI 
rule inhibits elastic scattering of the (j) on nucleons and 
non-strange mesons and resonances. Elastic-scattering 
cross-sections are dominated by Pomeron exchange and 
are expected to be on the order of 1 mb (see pj and 
references therein). In particular, from studies of pho- 
toproduction, total cross sections for the interaction of 
(j)'s with nucleons are estimated to be 8-12 mb over the 
range of kinetic energies from a few hundred MeV up to 
several GeV, and dominated by OZI-allowed absorptive 
processes — elastic scattering of (j)'s on nucleons accounts 
for about 1 mb of this cross section . It might there- 
fore be expected that (/)'s decouple rather early from the 
nuclear fireball, which is composed of nucleon resonances, 
pions, and non-strange meson resonances. If this were 
the dominant effect, the (j) inverse slope (T) would reflect 
the original "temperature" when the (/)'s are produced, 
without the enhancement from collective transverse flow 
developed in the late hadronic stages j4|]. In this case, 
T for the ip's would be lower than that observed for other 
hadrons of similar mass (for example, p's). 

On the other hand, when the (j) is detected via its decay 
into kaons, the rescattering of the daughter kaons from 
(/)'s that decay inside the nuclear fireball would generate 
a relative depletion in the (j) yield at low pt- This expla- 
nation has been proposed to explain the difference 
in the values of T for (p mesons obtained by experiments 
NA49 and NA50, which observed the (f> via its decay into 
K^K~ and iJ.'^fi^ final states, respectively. Given the 
short lifetime of the (cr = 45 fm/c), one might expect 
that this effect can give rise to a significant depletion in 
the mt spectrum only at low rrit- To illustrate, we have 
performed a simple simulation in which (jj's are generated 
at the center of a "fireball" according to an exponen- 
tial rrit distribution with T = 250 MeV/c^, and any 0's 
that decay within 7 fm of the origin are suppressed. For 
rrit — iTT-o > 0.5 GeV/c^ the resulting depression of the mt 



spectrum is nearly uniform; a deficit in the spectrum of 
more than 10% is noted only for TO( — mg < 0.25 GeV/c^, 
i.e., only at the edge of the E917 acceptance. (A calcu- 
lation performed with RQMD for Pb -I- Pb collisions at 
the SPS gives similar results [^.) Conceivably, however, 
in-medium effects could broaden the (f> enough to cause 
a more significant fraction of 0's to decay inside the fire- 
ball, leading to an increase in the apparent value of T in 
fits to the mt spectra (see Ref. |^| for an example with 
reference to the NA49 and NA50 data). 

For the comparison of the transverse-mass spectra for 
various particle species, we choose to use the mean trans- 
verse mass, (mt). This is because a variety of forms are 
used to fit and parameterize the transverse-mass spec- 
tra in the existing measurements of tt"*" 
p [|0[ |4|], d [|7), and A production at midrapidity in 
central Au -I- Au collisions at the AGS, including Boltz- 
mann, Boltzmann plus exponential, and mt-scaled expo- 
nential forms, in addition to the simple single-exponential 
form used to fit the mt spectra for 0's in this work. In all 
but two of the reports referenced above, (mt) is derived 
from the slope parameters for the chosen parameteriza- 
tion of the mt spectra, and quoted in place of the slope 
parameters themselves. The exceptions are the and 
A data. In Ref. the mt spectra for kaons are fit with 
exponential distributions, and values of T are reported. 
We convert to (mt) using the expression 



(mt) 



mlc^ + 2moc^T + 2T'^ 



m^c^ 



T 



(8) 



We use the same expression to obtain {mt) from the value 
of T obtained from the fit to the mt spectrum for the cf) 
presented in this work. In Ref. [|2|, the mt spectra for 
A's are fit with Boltzmann distributions, 



1 



dN/dy 



2Trc^mt dmtdy 



27r(rBmgc4 + 2r2TOoc2 + 2T|) 
(mt - mo)c^ 



xmtc exp 



Tb 



(9) 



and values for Tb are quoted. We convert to (mt) using 
the expression 



2moc^TB 



2Ti 



(10) 



The values of {mt — mo) for (f> mesons, together with 
the corresponding values for tt+'s, if+'s, p's, d's, and A's 
are plotted as a function of species mass in Fig. |. We 
note that the centrality selection differs somewhat from 
species to species in the available data: top 3% for 7r+, 
p (E866), and_ d; top 5% for p (E917), (j) and K+; and 
top 12% for A. The rapidity coverage varies similarly: 
for the 7r+ data it is 1.3 < ?/ < 1.4; for the and 
(j) data it is 1.4 < ?/ < 1.6, for the p and d data it is 
1.4 < ?; < 1.5, and for the A data it is 1.0 < y < 1.4. 
(For reference, y^N = 1-61 for 11.7^GeV/c Au + Au 
collisions.) The dashed line in Fig. |^ shows a fit to 
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FIG. 7: The mean transverse mass {{rrit — mo)) for various 
particle species as a function of the species mass at midra- 
pidity for central Au + Au collisions at 11.7 GeV/c per nu- 
cleon as measured by E866 and E917. The dashed line is a 
fit assuming exponential transverse-mass spectra and a lin- 
ear relationship between inverse slope (T) and particle mass 
(mo), as described in the text. Values of T for particles other 
than the are taken from Refs. |, |l2[ |30[ H, ^ . The pro- 
ton data point from E866 (E917) is presented with a shift of 
-30 MeV/c^ (-1-30 MeV/c^) in mass for clarity. 



the values for {mt — mo) using the form of Eq. (^) with 
T = a + brriQC^ , so that a and b are the free parameters of 
the fit. Assuming that it is qualitatively valid to describe 
the midrapidity toj spectra for the various species by ex- 
ponential inverse-slope parameters, this fit illustrates the 
trend expected for a linear relationship between inverse 
slope (T) and particle mass (toq), as motivated by hy- 
drodynamic models. 

The value of {mt — mo) for the (f) falls slightly below 
the systematic trend observed for the other particles, but 
less so than for A's. There may be some suggestion of an 
effect from early freeze-out of the (j) and late development 
of at least some of the transverse flow, as described in 
Ref. ji^. However, the significance of this observation is 
limited by the precision of the measurement. 



B. Centrality dependence of <?!) production; 
comparison to other species 

As mentioned in Sec. |, processes such as pp — > ppip are 
suppressed by large threshold energies and the OZI effect. 
Thus a naive expectation assuming ordinary hadronic in- 
teractions is that if an enhancement in <f) production were 
observed in heavy- ion collisions at AGS energies, this en- 



hancement would result from secondary collisions, e.g. 
via channels such as K'^ K (pp or K^K~ (pp [g4|. If 
this were in fact the case, increasing KY and KK combi- 
natorics would bring about an increase of the (p/ K ratio 
for central collisions. 

On the other hand, the proposed "re-arrangement" 
or "shake-out" of an intrinsic ss component of the nu- 
cleon wave function in the non-perturbativc regime pro- 
vides a mechanism for (p production in N + N interac- 
tions that bypasses the effects of OZI suppression. If the 
ss component of the nucleon wave function were neg- 
atively polarized, the fact that, at threshold, the reac- 
tion pp pp(p must proceed with the pp initially in the 
■^5*1 spin state would imply that this mechanism should 
be particularly important at near-threshold energies (see 
Ref. Q and references therein, especially Q). Indeed, 
the DISTO collaboration has observed that in p + p col- 
lisions at ^/s = 2.90 GeV, only 83 MeV above threshold, 
the (p/K~ ratio is about unity, such that <p production 
represents an important contribution to the K~ yield at 
these energies Q. 

Along these lines, a possible mechanism has been pro- 
posed to explain overall strangeness production in N + A 
collisions within the framework of the additive quark 
model Q. According to this proposal, strange particles 
are born as strange-quark pairs from binary collisions 
of the projectile and target nucleons, with a probabil- 
ity proportional to the number of interacting constituent 
projectile quarks In p + Au collisions at 17.5 GeV/c, 
E910 has shown that the production of A's and Ks^s in- 
creases with the estimated number of binary collisions, 
V, suffered by the incident proton jl^. For v < "i, the 
increase in the yields is faster than expected from scaling 
oi p + p data by the number of total participants (i.e., 
Np+A/Np+p = ^{1 + v)), although bounded from above 
by linear scaling with v (i.e., Np^^/Np^p = v). 

If the same mechanism were responsible for the ob- 
served increase in <p production with centrality in A -t- A 
collisions, we would expect to observe similar scaling be- 
havior in our data. Specifically, we would expect approx- 
imately constant numbers of hadrons bearing strange 
quarks to be produced per binary collision. In the follow- 
ing, we therefore compare the yield of (/)'s to the yields of 
pions and kaons as a function of centrality. The yields of 
pions and kaons were obtained by the E866 collaboration 
from Au -I- Au collisions at the same beam energy used 
for the present measurement 

As seen from the left panel of Fig. ^, the (p/n ra- 
tio of fiducial yields shows a rise toward central colli- 
sions, which signals an enhancement in the production 
of (p mesons relative to that of non-strange tt mesons 
in central collisions. This enhancement is also clearly 
suggested by Fig. |[ since pion production is known to 
mainly come from resonance decay and secondary rescat- 
tering and scales linearly with (Npp). Due to the fact that 
the plotted points represent ratios of fiducial yields, some 
care should be taken when interpreting the exact form of 
the dependence. We note that this increase in the (p/n 
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FIG. 8: The (^/tt and (f)/ ratios as a function of (Npp). For 
the (f>, the fiducial yields from E917 in the rapidity interval 
1.2 < y < 1.6 are used. For tt's, the data have been taken 
from Ref. and represent fiducial yields in the rapidity 

interval 1.2 < y < 1.4. The value used for the tt yield is 1.125 
times the yield for 7r+, which is essentially the average yield 
for tt""" and tt" . For 's, the data have been taken from 
Ref. and represent fiducial yields in the rapidity interval 
1.2 < y < 1.6. 
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FIG. 9: The normalized fiducial yield of 4> mesons per num- 
ber of binary collisions, (A'coii), as a function of (A'coii). The 
accepted rapidity region is 1.2 < y < 1.6. 



ratio is qualitatively similar to the increase in the K/tt 
ratio with centrality observed by E866 in Au + Au colli- 
sions at the AGS Q. 

We next compare the degree of enhancement in the 
yields of 0's and K's, which both contain s quarks. We 
use the yield of to characterize the overall kaon pro- 
duction because the E866 data indicate that the / 
ratio is about 0.15 and independent of centrality. In the 
right plot of Fig. ||, the ratio of fiducial yields (j)/K'^ is 
plotted as a function of centrality, and shows no substan- 
tial variation. This implies that both and (or K^) 
possess a similar degree of enhanced production toward 
central collisions at AGS energies. As estimated from 
linear fits to the 0/_ftr+ ratios as a function of {Npp), an 
increase of up to 50% for the ratio in the most central 
bin cannot be definitively ruled out within a la range of 
the fit errors. However, any centrality dependence of the 
4>/K ratios appears to be weak. 

To look for scaling behavior similar to that observed 
in p + Au collisions by E910, we examine the dependence 
of the (f) yields in our data on the mean number of binary 
collisions, (A'^cou)- (Estimation of (iVcou) is discussed in 
Sec. |l[ the values of (iVcou) for each centrality bin are 
given in Table |.) The fiducial (j) yield normalized to 
(A'coii) is plotted versus (A'coii) in Fig. ||. The dependence 
of the fiducial (p yield on (A'coii) is consistent with scaling 
proportional to (A'coii). 

While the rapidity coverage of our measurement does 
not allow us to make a precise statement about the ab- 
solute value of the yield, it is possible to make an in- 
formed guess about the width of the </> rapidity distribu- 
tion. Both E859 g and NA49 ||, || have observed that 
the Gaussian width of the </> rapidity distribution is very 
similar to those of the and K~ rapidity distributions. 
In central Au -|- Au collisions at the AGS, the Gaussian 
widths of the rapidity distributions for and K~ are 



^--^0.9 and ^^0.7, respectively y, |25|. The narrow 
width has been attributed to the higher threshold and 
more restricted phase space for K~ production. Since 
the thresholds for K~ and (f> production in p + p col- 
lisions are very similar, we assume that the (jj and 
have Gaussian rapidity distributions with approximately 
equal widths in central Au -I- Au collisions. We then esti- 
mate that our fiducial yield corresponds to about 20% of 
the total yield. Using the parameterization described in 
Appendix ^ for the (f) yield in p + p collisions as a func- 
tion of center of mass energy, ^/s, the yield of (/)'s per 
binary collision in Au -I- Au collisions is about 50% of the 
(j) yield in p p collisions at this energy. This observa- 
tion seems contrary to the usual expectation of enhanced 
strangeness production in Au -I - Au collisions. We discuss 
this point further in Sec. IV C . 

Using E866 data on the total yields of and AT^ 
normalized to (Acoii), a similar scaling proportional to 
(A^cou) is observed. It appears as if the mechanisms for 
the production of AT+'s, Ar~'s, and 0's all have a similar 
dependence on the centrality of the collisions, and this 
dependence is consistent with a scaling with the number 
of binary collisions, (A'coii)- This observation suggests 
that hard binary collisions might play an important role 
in strangeness production in heavy- ion collisions. 

Our observation that (f> and kaon production scale sim- 
ilarly and faster than linearly with (Npp) therefore pro- 
vides an essential test of the details contained within 
rescattering models. The rescattering model that sug- 
gests the importance of contributions from KY and KK 
collisions in the increase of (f> production in central A -I- A 
collisions can be ruled out; if such contributions were 
dominant, combinatorial considerations would lead to the 
expectation that (j) production should increase faster than 
kaon production with centrality. This is incompatible 
with our observation that (p arid kaon production show 
similar scaling with centrality. 
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FIG. 10; Ratios of 0/(7r) in central heavy-ion and p + p in- 
teractions as a function of y^sjvjv (top panel) and double ra- 
tios of A -I- A over p + p (bottom panel). For the E917 and 
STAR points, fiducial yield ratios are plotted; for the other 
two points, total yield ratios are shown instead. The hatched 
area represents the ratio of total yields in p+p collisions, based 
on the parameterization discussed in Appendix The cross 
bars in the bottom panel indicate the contribution to the to- 
tal errors from the uncertainty on (f)/ (ji) in central heavy-ion 
interactions. 



C. Dependence of ^/tt, (jj/K^ , and 4>/K ratios on 
y/SNN in A + A and p + p reactions 

In order to further explore the mechanism responsible 
for the approximate scaling of the yield with {Ncow) 
observed in our data, we compare our results with other 
measurements of 4> production in heavy-ion collisions at 
different reaction energies. 

The excitation functions of the ^/tt, (j)/K'^, and 
(j)/K~ ratios in central heavy-ion collisions are shown 
in Figs. |l^, |ll|, and ^ The four points correspond 
to Au + Au collisions at the AGS (this measurement. 



FIG. 11: Same as Fig. for the ratios of 



= 4.87 GeV) i, p 
{^/I^ = 5.39 GeV) [In 




Au collisions at the AGS 
Pb + Pb collisions at the 
3| , and Au + Au collisions 
|, These plots must 
Most obviously, the coUi- 



SPS {y/IWN = 17.27 GeV 
at RHIC {y/sWN = 130 GeV) 
be interpreted with some care, 
sional system is different in each case; in particular, the 
various particle yields are not guaranteed to scale in the 
same way when passing from the Si + Au system to the 
heavier systems. In addition, the points for Au-fAu colli- 
sions at the AGS and RHIC represent fiducial-yield ratios 
(the AGS point is for the fiducial yield over 1.2 < y < 1.6, 
and the RHIC point is for the central unit of rapidity), 
while the other two points are ratios of yields over all 
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phase space. Nor have we applied any corrections in the 
comparisons with p+p collisions to take into account the 
isospin averaging of yields from p + p, n + n, and p + n 
collisions ^ Nevertheless, two observations can 
be made. 

Our first observation is that the ^/tt ratios for A -I- A 
collisions are notably enhanced with respect to their val- 
ues for p + p collisions for all points except at the highest 
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FIG. 13: Same as Fig. |l^ for the ratios of (j>/{Ncoii), where 
(A'coii) is the number of binary collisions. For the E917 data 
point, an estimated total (p yield is plotted as described in the 
text. 



The ratio 4>/{^coi\) for central A + A collisions at the 
AGS and SPS is plotted in Fig. |l^. For the comparison, 
we are forced to assume a value for the 4> rapidity width 
for our measurement. A width of ay — 0.71 equal to the 
measured K~ width is used, as explained in the previ- 
ous section, and we plot our point with an additional 20% 
systematic uncertainty corresponding to a range of values 
for ay from 0.6 (overlap of and K~ rapidity distri- 
butions) to 0.9 {K+ rapidity distribution). The RHIC 
point is not included in this comparison, because there is 
no reasonable way to extrapolate the fiducial yield to all 
of phase space. The ratio (^/(iVcoii))A-i-A/(0/(A'"coii))p-i-p 
is consistent with ^0.5-0.7 for all three points from the 
AGS and SPS. Our comparison is by no means precise. 
However, it does seem that the yield of (/)'s per binary col- 
lision, modulo the effects of threshold and center-of-mass 
energy dependence, is approximately constant across the 
three collisional systems studied. Furthermore, instead 
of there being an enhancement in the yield per binary 
collision in heavy- ion collisions, this double ratio is less 
than 1. This might reflect the effect of (j) absorption in 
the nuclear fireball, as the inelastic cross section of the 
(j) on nucleons (~8~10 mb) is a significant component of 
the total interaction cross section (^--^8-12 mb) H, |55[|. 



SUMMARY 



energy in Fig. |l0|. In p + p collisions, (^/ir increases with 
y/^NN, in large part because of the larger production 
threshold for (p mesons. In heavy-ion collisions, the (I)/tt 
ratios are enhanced by factors of 2-3, at least up to SPS 
energies, but seem to show the same energy dependence 
observed in p + p collisions, at least up to SPS energies. 
Such an enhancement in the central A -I- A collisions could 
be interpreted as due to the different scaling behaviors 
of (j) and TT with centrality. As discussed in Sec. IV B , 
the yield of iji's scales with (iVcou) while that of tt's scales 
with (Npp). The ratio of (A^cou) to (A^pp) is around 1 
in peripheral collisions, which are similar to p + p col- 
lisions, and becomes larger than 1 in central collisions. 
At RHIC, the ^/tt ratio for central Au + Au collisions 
is about the same as that from the parameterization of 
the p + p data. It might be speculated that the primary 
production channels for (p's and/or tt's at RHIC energies 
are different from those at AGS and SPS energies. 

The second observation is that the 4>/K ratios for A-t-A 
collisions are only marginally enhanced with respect to 
their values for p + p collisions, and show relatively little 
variation as a function of energy as seen in Fig. |llj and 
Fig. ^ Our parameterization of the p + p cross section 
for production suggests that the (j)/K~ ratio increases 
as ^snn approaches the threshold value from above for 
the reaction pp — > pp(j). The (j)/ K~ ratios in heavy-ion 
data seem to follow the energy dependence observed in 
p + p collisions rather reasonably. 

It is interesting to see if the data points from other sys- 
tems and energies obey the scaling with (A^coii) as well. 



In conclusion, we have studied (j) production in Au-|-Au 
collisions at 11.7AGeV/c around midrapidity as a func- 
tion of collision centrality. The yield per projectile partic- 
ipant shows a steady rise toward central collisions. This 
enhanced production in central collisions is stronger than 
that of non-strange tt mesons as seen from the increasing 
(f)/TT ratio with centrality. The ratios (j)/K~^ and (t>/K~ 
are approximately constant with (iVpp). The yield of 0's, 
like the yields of and K~ ,is seen to scale with ( A^cou) , 
the number of binary collisions, and this observation is in- 
compatible with predictions from some rescattering mod- 
els of (j) production, in which (p production increases faster 
than kaon production with centrality due to combinato- 
rial considerations. Finally the yield of (j)'s per binary 
collision in A + A collisions is about 50-70% of the 
yield in p -f p collisions at AGS and SPS energies. That 
the yield of (/)'s per binary collision in A + A collisions 
is smaller than that in p + p collisions might signal the 
effect of (j) absorption by nucleons in heavy-ion collisions. 
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APPENDIX A: PARAMETERIZATION OF (/> 
PRODUCTION AS A FUNCTION OF IN p + p 
COLLISIONS 

The (j) yield in p + p colhsions provides a useful bench- 
mark for the interpretation of (j) yields in heavy-ion col- 



lisions. There are various existing measurements of the 
inclusive total cross section for (j) production in p + p 
collisions for 5 < y/s < 60 GeV ||, Q that can be 
extrapolated to the values of ^/snn for which heavy-ion 
data exist. (Note that the point at -^/s = 2.90 GcV from 
Ref. p8| , while technically an exclusive measurement in 
the pp — > ppK~^K~ channel, is also an inclusive measure- 
ment, since no other channels with a or in the final 
state are kinematically allowed at this energy.) We have 
fit these measurements with a form used to parameter- 
ize results on vector-meson production cross sections in 
p + p collisions obtained from one-pion exchange calcula- 
tions Q and the Lund String Model ||§ : 



aipp ^ <j)X) = a(l - xfx'' 

^ = -^thresh/ ^ 

sthrcsh = 8.38 GeV^ (Al) 



The results of our fit are shown in Fig. |l^. We obtain 
a = TAtll Mb, b = 2.05t!^:i8, and c = 1.56l|]:^?. For 
comparison with heavy-ion data, we obtain yields by di- 
viding the (f) cross sections for each value of ^/snn by the 
inelastic p + p cross section. For Au -I- Au collisions at 
the AGS, y'SNN — 4.87 GeV, and our parameterization 
of the p -hp data gives a yield of 0.00169 ±0.00042 at this 
energy. 

For K and tt production in p + p reactions, we use the 
multiplicity parameterizations in Ref. . The smallest 
■\/s for the ranges of parameterization for tt, if + and K~ 
are 3.0, 2.98, and 5.03 GeV, respectively. 
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